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The relativistic form of the Vlasov equation is used to solve the problem of transmission and 
reflection of normally incident electromagnetic waves. Transmission and reflection coefficients 
are derived for a plasma half-space and for a plasma slab, using the assumption that electrons are 
specularly reflected at the boundaries. These coefficients are functions of the ratio of the electron 
thermal speed to the vacuum speed of light, indicating that temperature dependence in these cases 
is a relativistic effect. Nevertheless, it is seen that the nonrelativistic limits of these coefficients 
differ from those obtained by using the cold-plasma equations since the transition to the zero-temper- 
ature limit has a nonuniform character. 



1. Introduction 

Investigations of transverse electromagnetic wave 
propagation in isotropic plasmas are usually based on 
the cold-plasma equations. Since such an approach 
reveals that the phase velocity of the waves is larger 
than the vacuum speed of light c, and therefore much 
larger than physically admissible thermal velocities 
of the plasma particles, it is concluded a posteriori 
that thermal effects are negligible. The validity of 
the cold-plasma approximation for the case of a trans- 
verse wave propagating in an unbounded collisionless 
plasma has been demonstrated by Buneman, [1961a], 
who derived a relativistically invariant dispersion 
relation by solving the equations of motion and con- 
servation of the different electron streams that are 
superimposed on each other in their combined fields. 
In contrast to the case of longitudinal waves, in which 
the inclusion of temperature effects is necessary to 
obtain a dispersion relation, Buneman's dispersion 
relation for transverse waves contains a thermal cor- 
rection to the cold plasma result that is negligible for 
laboratory temperatures. 

The large phase-velocity argument which motivates 
the use of cold-plasma formulas applies only to regions 
of the plasma which are far away from boundaries. 
In investigations of longitudinal wave problems, a 
kinetic approach has revealed the existence of a short 
range disturbance which is a superposition of modes 
which have phase velocities ranging from zero to the 
vacuum speed of light [Taylor, 1963]. It is expedient, 
then, to use a relativistic kinetic analysis to investigate 
the effect of boundaries on the transmission of trans- 
verse waves, since the slow components of the field 
provide a mechanism for the exchange of energy be- 
tween the field and the plasma even in the absence of 
collisions. In this paper, the approach previously 
developed for longitudinal waves is extended to the 
transverse wave problem. We will obtain a dispersion 
relation for the propagating component of the field 
which is identical to that of Buneman, and also a 



tractable expression for the short range disturbance 
produced by the effect of the boundary on the elec- 
trons. We will then define transmission and reflec- 
tion coefficients for the plasma half-space, and contrast 
these with the results obtained with the cold plasma 
approximation. Finally, we will treat the problem of 
electromagnetic wave propagation in a plasma slab. 
In both of these cases temperature dependence of the 
coefficients vanishes when the electron thermal veloc- 
ity may be neglected when compared to the speed of 
light, but the limiting results differ from those obtained 
by using the cold-plasma equations. 

It should be noted that as in most treatments of 
problems of this type, we will assume specular reflec- 
tion of electrons at the boundary. This assumption is 
based on the existence of a potential sheath [Bohm 
and Gross, 1950]. 

2. Formulation of the Problem 

We erect a rectangular coordinate system so that 
the plasma is contained in the half-space z > 0. An 
incident transverse wave with frequency cu propagates 
in the vacuum in the positive z-direction so that its 
electric field is alined with the positive ^-direction. 
We will determine the fields in the plasma and the 
amplitude of a wave reflected into the vacuum. Be- 
cause we wish to solve this problem by means of a 
Fourier transformation, we use a mathematical arti- 
fice which was introduced by Shafranov [1958]. We 
consider an infinite plasma with a sheet current source 
at the origin z = 0. Setting E(z) = e x E(z)e~ i<l>t , we have 



dz 2 



E + — E+ iajfJLtf = — A8{z) 



(1) 



where j is the plasma current and A is an arbitrary 
constant. The solution of (1) for all z is certainly a 
permissible solution for the half-space z > of our 
original problem. We need only require that at z = 0, 
it connects a physically proper way to the vacuum field 
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in the left half-space, which means that the constant 
A must be chosen so that the electric and magnetic 
fields are continuous across the boundary. 
The plasma current is given by 



■/ 



j = e d 3 u/3cujf, 



(2) 



where / is the solution of the relativistic collisionless 
kinetic equation [Clemmow and Willson, 1956] , 



. iwf+ p CUz f r -^ E ^F =l+u >. 



(3) 



The notation here is the same as that used previously 
in the longitudinal wave problem [Taylor, 1963]. 
Proceeding in a straightforward manner analogous to 
that of the longitudinal wave treatment, and assuming 
specular reflection of electrons at the boundary, we 
obtain 



dz 2 ' 



-^)e- j^ Kj(\z-z'\)E(z f )dz' 



JO 



Ki(z+z')E(z')dz'=-A8(z), (4) 



written as a sum of two components E\ and E-i, where 



E,=-iA 2 Res 



((o/cf-k*-Kj{k)} k = k . 



E->(z) = + : 



I 



e ikz dk 



2Tr)c(a>lc) 2 -k 2 -K J {k) 



(9) 



(10) 



Each kj in the summation of (9) is an isolated pole of 
the integrand of (7), and the path C proceeds clock- 
wise around the branch line in the upper half of the 
Zc-plane. 

3. Dispersion Relation 

To calculate E\ we must find the roots of the dis- 
persion relation 



~k*-Kj(k) = 0. 



We set n = kclo) to obtain 



l-n 2 -^/(/i) = 0, 
or 



(11) 



(12) 



where 



JW*)-= 



r+oo r+oc fx u df(\ 

dii x du y diiz — —^- 

J-x J-oo ' Jo u z 6u x 



exp tej- 



(5) 



Now we can consistently set E(z) = E(—z) to obtain 

{^+^) E -J m K,(\z-z'\)E(z'W=-A8(z). (6) 

A Fourier transformation of (6) yields 
„, , A [^ e ikz dk 



-KAk) 



where 



K T (k) = ^ [ cPu 



fSu x 



Ms. 

du x 



kc 



Pu z -l 



(8) 



We have set co= \co\e i8 so that the integrand is defined 

everywhere in the /c-plane except at those points on 

the branch line k= \k\e i8 for which \k\ > I — I . 

c 

We may simplify the integral in (7) by deforming the 

contour of integration in the A>plane. The result of 

such a procedure is that the electric field may be 



^here 



Kn)-- 



-\d>u 



nfiiiz — 1 



It is shown in the appendix that this relation may be 
written to first order in the ratio of electron thermal 
energy to electron rest energy as 



1 



1+( * 2 - 5/2) ii]= ' 



(13) 



which is the same as that obtained for the unbounded 
plasma by Buneman [1961a], who used a different 
approach. If we give a> p the same small phase angle 
d that we have given a>, then (13) is a linear equation 
in n 2 with real coefficients, and it has one real root. 
It can be shown that no roots have been lost in the 
approximation by examining the exact dispersion 
relation given by (12). This means that in the evalua- 
tion of Ei we need consider only one pole in the k- 
plane, and that this pole lies on the positive real or 
the positive imaginary axis in the collisionless limit. 
A simple computation yields 



Ei(z)- 



iAe lk v z 

Cittn 



(14) 



where 



k p = co/c 



/l-(a>> 2 )(l-|/Cr/mc 2 ) 
l+(a>l/<o>)(KTIm<*) ' 



(15) 
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It follows from (15) that whenever co 2 



<«$ 



,1-1 



KT/mc 2 ) we obtain a nonpropagating exponentially 

damped disturbance for E\, while for co>co 2 (1— ■ — 

KT/mc 2 ) we obtain a propagating unattenuated wave. 
This conclusion represents a relativistic correction to 
the well-known cold plasma result. Maxwell's 
equations lead to a magnetic field B\ associated with 
the electric field that we have computed, and it may 
be written 



Bx 



iyBi e y 2w e" 



(16) 



4. Short Range Component of the Plasma 
Field 

We may evaluate the field E> by the extension of 
Landau's method [Landau, 1946] which was previously 
developed for the relativistic treatment of longitudinal 
waves [Taylor, 1963]. We allow 8 to go to zero so 
that the branch cut swings onto the real axis. The 
integral around the cut on the real axis may be re- 
duced to 



E2(Z)=~~ 



A 
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/: 

Jtolc 



[K-(k)-K + (k)]e ik *dk. 



[ (cole) 2 - k 2 - K (k) ] [ (co/c) 2 -k 2 - K + (k) ] ' 



(17) 



from those obtained by the cold plasma approximation 
even for very small temperatures. This is due to the 
behavior of the short-range component of the electric 
field — its amplitude and its range are both proportional 
to (KT/mc 2 ) 112 , so that the magnetic field which is 
induced by it has a finite limit as T— * 0. 



5. Transmission and Reflection Coefficients 
for a Half-Space 

Combining the results of sections 3 and 4, we find 
that to zero order in v ih /c the fields at the wall are 
given by 



*«»-£. *<o)-j£(i+«), 



(21) 



wnere a~- 



V2< 



The parameter a represents the 



correction to the results obtained by using the cold- 
plasma equations — the latter results are given by (21) 
when a is set equal to zero. If we define the trans- 
mission and reflection coefficients T and R by the 
equations 



E T (z)=TEie ik p z e z , 
E R (z)=RE I e- ik p*e z , 



(22) 



where 



CO~CO f+oc 

K±{k) = -R-\ dux 

C 2 J -oo 



\Z dUv \c ± 



du z 

Pu x (df {) lduj-) 
kc f3u z — co 



(18) 



The path C+ is a semicircle in the upper half of the 
a^-plane, while the path C. is a semicircle in the lower 
half of the a^-plane. It may be shown, using the 
method of steepest descents that £2(2) and B 2 {z) have 
negligible amplitudes several Debye lengths away from 
the boundary, which means that the range of these 
components is proportional to the square root of the 
temperature of the plasma. The amplitudes of these 
fields at the interface may be computed by setting 
z = in the integrands of their representations and 
expanding in powers of (KT/mc 2 ) . We obtain 



£ , 2 (0)=0[(/^r/mc 2 ) 1 / 2 ] 



(19) 



fi 2 (0) = 



iAt 



2V2o 



- + 0[(AT/mc 2 )]. 



These results, combined with (14), (15), and (16) 
indicate that there is a nonuniform transition to the 
zero-temperature case, so that the transmission and 
reflection coefficients of the plasma interface differ 



where Ei is the amplitude of the incident electric field, 
Er(z) is the electric field in the plasma far from the 
boundary, and Er(z) is the field reflected into the 
vacuum, we obtain 



1+y ' 



(23) 



R=- 



/here 



y = Vl-(co p /co) 2 (1 + a). 



If we define P/, Pt, and Pr as the incident, trans- 
mitted and reflected power in the electromagnetic 
field, we obtain 




co > co u 



co < co t) 



(24) 



(20) and 




CO < C0 D 



CO > COp) 
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The cold plas ma results are obtained by setting 

y=Jl- (^f in (22), (23), and (24). It is interesting 

to note that the sum of the transmitted power and the 
reflected power is not equal to the incident power 
when a) > co p , indicating that some energy is delivered 
to the electrons in the plasma. The fraction of inci- 
dent power delivered to the plasma may easily be 
shown to be (4ya) (1 -hy)~ 2 (1 + a) -1 . (The concept 
of kinetic power has been discussed by Buneman 
[1961b] and by authors listed therein.) 



6. Transmission and Reflection Coefficients 
for a Plasma Slab 

In sections 3 and 4 we found a solution to the one- 
boundary problem in which there is no wave incident 
on the boundary from within the plasma. We j?an 
generalize this solution by adding the solution E of 
the equations 



. r , df e fr df 

— icx>f+ficu z — = -E- — , 

dz mc ou x 



(25) 



and 



dz 2 



E + (o 2 lc 2 E + io)fjLQj = 0, 



(26) 



which is consistent with the requirement of specular 
reflection. It is easily shown that such a solution 
is given by an electric field which is proportional to 
cos k p z, where k p is defined by (15). Several Debye 
lengths from the wall, where the short range component 
of the fields can be neglected, the general fields are 
given by 



E(z) = (A + BI2)e ik P z + B/2e- ik P z 
B(z) = k p /o) [(A + B/2)eV - B/2c" V] , 



(27) 



where A and B are arbitrary constants. At the 
boundary, where the short range component must 
be included, the fields are given by 



E(0)=A + B 
B(0) = kplmA(l + a). 
If we write the vacuum field as 

Eie ik o z + E R e- ik o z 

and the field inside the plasma as 

E(z) = E PI e~ ik p z + E PR e +ik P z , 



(28) 



the boundary conditions become 

*.=*=!- • 2 



■ Epi + 
1 + 7 1 + 7 



E,, 



Er — — — Epi - 
1 + 7 



1 + 7 



Er 



(29) 



If we use these boundary conditions to solve the prob- 
lem of transmission by a plasma bounded by two 
specularly reflecting walls of separation d, and then 
define E(z) = Eie ik & + RE ie - ik o z , z<0 and E(z) 
= TEie ik o z , z > d, we find 



T= 



47 



(1 + 7) 2 e~ ik P d - (1 - 7) 2 e ik P d 



R = 



1 + 7 



(1 - Te ik P d ). 



(30) 



The ordinary cold-plasma results are again obtained by 
setting a = 0, which leads to y = Vl — a) 2 /(o 2 . It can 
easily be verified that |T| 2 + \R\ 2 = 1, indicating that 
there is no net exchange of energy between the field 
and the electrons in the zero temperature limit, 
but that the zero-temperature limit of transmission 
and reflection coefficients for the slab obtained from 
kinetic theory differ from those obtained from the cold 
plasma formalism due to the presence of a in the 
definition of 7. 



7. Discussion 

We have explored the consequences of the relativ- 
istic form of the Vlasov equation when specular 
reflection of electrons at boundaries is assumed. 
This problem has been formulated for a half-space 
by Felderhof [1963], who used a nonrelativistic normal 
mode analysis of the type introduced by Van Kampen 
for the study of longitudinal waves. Although 
Felderhof did not obtain results in closed form, he 
was able to conclude that the electromagnetic field 
can exchange energy with the plasma electrons, and 
that there is an unattenuated propagating wave in the 
plasma at sufficiently high frequencies, as we have seen 
in section 5. When the problem is solved for the 
plasma slab, however, we have seen that the net 
energy exchange between the field and the electrons is 
zero as long as the plasma is rigorously collisionless. 

It is interesting to note that the actual temperature 
dependence of our coefficients is of order (KTImc 2 ) 112 , 
so that temperature effects are usually negligible, but 
that the transition to the zero temperature limit is 
nonuniform. 

It can easily be shown that the coefficients of trans- 
mission and reflection for a plasma in a uniform 
magnetic field are the same as those we have derived 
if the electric field vector of the incident wave has the 
same orientation as the fixed magnetic field. 
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8. Appendix 

To obtain the approximate eq (13) from the exact 
eq (12), we must evaluate the expression I(n). The 
technique given here is a generalization of the work of 
Clemmow and Willson [1956] for the case of longi- 
tudinal waves. Transforming to spherical coordinates, 
we have 



sin 3 6 cos 2 <f) df 

cos 0-1 ^du" ( A1 > 



Jo Jo Jo n pu 

I 

Now we make the transformation f = cos 6 to obtain 

(A2) 

'"- J -1 "S ~ LIP" 

We define the functions 



)o du J -i nlj—lipu 



2h(u) = 



d€ 



and obtain 



u) J_ in £-l/jBa 

d_ht \-. ~P 
du w 1 + uHl-n 2 ) 



du &U) l + u 2 (\-n 2 ) 



(A4) 



Nc 



defir 



d u ., dfo(u) 

-j- H(u) = u l —j — 
du du 



-fG(u) = (l + u>)4^- 
du du 



(A5) 



that 



/=2tt f " du [( j^#(")) h(u)+ (^G(u)W) 



(A6) 



Integration by parts leads to 

1-9 (*J (l + u 2 )H(u)-u 2 G(u) 

7-277 J„ d "(l + « 2 ) 3/2 [l + (l-« 2 )" 2 ] lA7) 

Now we take the electron distribution function to be 



fo- 



277/ 



exp 



Xm 2 \ = mc 2 

2 /' xr 



(A8) 



Expanding the denominator of the integrand and re- 
taining the first two terms we have 



r =27T f* 

Jo 



du 



(A9) 



H(u)[l + (n 2 -3/2)u 2 ] 

-u 2 G(u)[\+(n 2 -5l2)u 2 ] 
Partial integration and use of (A5) lead to 

/ = 4t7 f da/o(«) [« 2 + (j- 5/6) u 4 ] ■ (A 10) 
Evaluating the integrals we have 



mr> z 



(All) 



(A3) and insertion of this expression into (12) yields eq (13). 
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